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ABSTRACT 

High resolution N-body simulations using different codes and initial condition techniques re- 
veal two different behaviours for the rotation frequency of transient spiral arms like structures. 
Whereas unbarred disks present spiral arms nearly corotating with disk particles, strong barred 
models (bulged or bulge-less) quickly develop a bar-spiral structure dominant in density, with 
a pattern speed almost constant in radius. As the bar strength decreases the arm departs from 
bar rigid rotation and behaves similar to the unbarred case. In strong barred models we detect 
in the frequency space other subdominant and slower modes at large radii, in agreement with 
previous studies, however we also detect them in the configuration space. We propose that 
the distinctive behaviour of the dominant spiral modes can be exploited in order to constraint 
the nature of Galactic spiral arms by the astrometric survey GAIA and by 2-D spectroscopic 
surveys like CALIFA and MANGA in external galaxies. 
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1 INTRODUCTION 

Since the early seventies, it has been suggested that the dynamics 
driven by bars and spirals have profound consequences on the kine- 
matic and structu r al evolution of galactic disks (e.g . Miller et al. 
' Hohllll97ll: lAthanassoulal [l980l : ISellwood & Athanassoula 
Friedli & Bend 1993 ). More recently, stellar radial migra- 
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tion in disk galaxies has been recognized as a critical component of 
disk galaxy evolution. This process may drastically alter our view 
of the connection between the present-day phase space and chem- 
ical distri butions of stars and the pro cesses of disk formation and 
evolution. ISellwood & Birmevl ( l2002h set up the dynamical frame- 
work of this process through the effects of transient spiral struc- 
ture, which seems to be a crucial process to underst and the solar 
neighb orhood observat ions as sugges t ed orig i nally bvlWielen et al 
' l996h . Authors such iRoskar et al.l (I2008I). ISchonrich & Binnev 



ISellwoo j I2OOOI: iD'Onghia et al.l [2013h. or the rnanifold theory 
llRomero-Gomez et al.l 20071 : iTsoutsis et ak I I2OO9I: lAthanassoulal 
I2OI2I) . 

In the quest of dynamical models not limited to quasilinear ap- 
proximations or steady state, N-body simulations have been used 
to understand the origin and ev olution of sp i ral structures. After 
pionee ring studies like Miller & Smithl ( [l97a) : ISellwood & Sparkd 
i ll988h : lRautiainen & Sale a99^ only recently appeared another 
boom of papers likely as a result of progress in the computational 
resources and codes. This re-ignition of simulations based studies 
opened-up an interesting debate about the pos sible corotating na- 



ture of spiral patter ns with disk particles (e.g. lOuillen et al 
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Baba etalj2013l) . As an example. lOrand et al.l ( l2012al lbl) computed 



20091) , iRoskar et all ( 1201 ll) and iMinchev et al.l ( l2012h . among 



others, revived the study of spiral arms and bars as triggers of 
stellar radial migration. Despite these numerous studies, funda- 
mental questions arise such as: what is the nature of the spi- 
rals? or, how is their pattern speed related to the motion of the 
stellar component? Several models have been proposed up to 
now, from the c lassical Tight-Winding Approximation (TWA, e.g. 
iBinnev & Tremaine 2008.) to th e mechanisms proposed to account 
for self-excited spiral patterns (lToomrelll990HBertin & Linll 19961 : 



collisionless N-body and smoothed particle hydrodynamics (SPH) 
simulations of disk galaxies and illus trated that transient spiral fea- 
tures appear to corotate with the disk. lComparetta & Ouillenl ( l2012h 
also suggested that some short lived features arising from con- 
structive interference between longer lived modes, i.e. fast bar and 
slowly moving spiral pattern modes, can be nearly corotating with 
the disk. This corotating nature wo uld have consequen ces on the 
stellar radial migration mechanisms jGrand et al.ll2012al) . 
From the observational point of view, the nature of spiral arms is 
also far from being clear (Sheth & Rossi 2010; Foyle et al. 20ll|; 
iFerreras et al.l I2OI2I) . Nowadays, only weak observational con- 
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straints are available to answer these questions. Some constraints 
come out from direct measurement o f the rotation frequ encies ra- 
dial variation in external galaxies (e.g. lMeidt et al .120091) , a method 
with a strong potential but currently applied only to a handful of 
galaxies, or from indirect m easurements as the one proposed by 
iMartmez-Garcia et al.l ( |2009|) . Other constraints to spiral arms' na- 
ture are based on small scale stellar kinematic substruc t ure an alysis 
in the Milky Way (MW) disk (e.g. lAntoia etallllOl ll , [201^ . Fur- 
thermore, although it is now commonly accepted that the MW is a 
barred galaxy, it is not clear how the arms and bar are related or if 
they are connected at all. The rotation frequency of both, bar and 
spiral arms, seem to show different values, none in mutual coro- 
tation or in corotatio n with galactic disk material , leavi ng unclear 
spiral arms' nature jMartos et aLll2004 lOerhardllioill) . It is fair 
to say that currently, gal actic arms' pattern speed estimations are 
mostly model dependent jMartos et al ]|2004l). The situation is also 
not clear for external galaxies ( Buta et all 2005h . The interpreta- 
tion of observational results is based on previous theoretical stud- 
ies and also on numerical simulations JSellwood & Sparkdligssl : 
iRautiainen & Sa lo 1999). Some of these studies adopted simplifi- 
cations like 2D N-body models and a rigid halo or no halo at all. It 
is not clear if such assumptio ns may affect t he gen erality of their 
conclusions as is suggested bv I Athanassoulal l l2002h for the case of 
bar growth. It is also not obvio us if the c orotating nature of spi- 
ral arms recently suggested bv lOrand et al . (2012aJ is valid also 
in models with different structure (i.e. models with strong/week 
bar and/or bulge). Furthermore, classical methods to estimate disk 
modes' rotation frequencies like time Fourier spectrograms applied 
to simulations, have been claimed to suffer from biases in models 
with multiple o r week spiral arms, hampering estimations of arms' 
pattern speeds JOrand et"al .l20I2al) . 

In this paper we analyze 3D galaxy models, with different stel- 
lar/dark structure, using live halos and with enough mass, force 
and time resolutions to accurately describe the internal disk kine- 
matics. We performed some testing on the results dependence 
on codes and initial conditions techniques. It should be empha- 
sized that the N-body simulations presented here are purely stel- 
lar. From the first attem pts to simulate gas in barred galaxies 
dSanders & Huntlevlll976l) . it has become clear that gas influences 
disk stellar dynamics, by even chang ing the live time of bars 
teird et al.ll20I2l : lDi Matteo et alj|20T3n . Gas makes disks to show 
more complicated and well-defined moiphologies. As an exam- 
ple, i nner rings, such the G alactic Molecular Ring observed in our 
MW jClemens et al.lll98^, h ave been the subject of numerous ob- 
servational teuta etal.ll2004l) and theoretica l investigations (e.g. 
iBvrd et al.l2006l : iRomero-Gomez et al l|). The influence of the 
gaseous component and its subgrid physics in our simulations is 
now un der investigation using the ART co de in the hydrodynamic 
version ( lKravtsovll2003l : IColi'n et alj201(]|) . 

The aim of the first study presented here is to revisit the possible 
correlation between spiral arm kinematics and their nature using 
the purely stellar component. In Section |2]we describe our N-body 
simulations, carried out using two well known N-body codes, ART 
and GADGETS. In Section [3] we describe the techniques used to 
derive the rotation frequencies, whereas the results and conclusions 
are presented in Section |4]and |5] respectively. 



2 MODELS AND SIMULATIONS 

We have performed coUisionless N-body simulations with ART 
( lKravtsov et"anil997h and GADGETS (last described in ISpringell 



Parameter 


Bl/5 


Gl 


Ul/5 


Disk mass (10^° Mq) 


5.0 


6.01 


3.75 


Halo mass (lO"*^^ Mq) 


1.38 


0.66 


1.5 


Disk exp. length (kpc) 


3.86 


3.0 


4.0 


Disk exp. height Z^i (kpc) 


0.2 


0.2 


0.2 


Halo NFW i?d (kpc) 


29.19 


14.4 


16.61 


Halo concentration 


10 


10.4 


18 


Halo DM species 


6/7 


1 


6/7 


N.d,sfe(+N.i,„;ge) {10^) 


1.0/5.0 


0.5(+0.5) 


1.0/5.0 


Neff (10^) 


2.86/13.8 


0.15 


4.1/20.0 


Min. time step (10* yr) 


3.2/1.6 


4.4 


7.9/3.1 


Spatial Resolution (pc) 


44.0/11.0 


35.0 


11.0/11.0 


Total integration t (Gyr) 


4.6/2.8 


1.4 


3.2/2.8 



Table 1. Parameters of the simulations: Column indicated by Bl/5 presents 
the parameters for models with 1 million (Bl) and 5 million (B5) star parti- 
cles in the disk, simulated using ART code, idem for unbaiTcd models Ul/5 
in the third column and for the barred model Gl, using GADGET3 code, in 
the second column. 



I2OO5I) codes. Tests of consistency of bo th codes, applied to 
dynam i cs of barred galaxi es can be found in Valenzuela & Klvpinl 
( l2003h . lColin et all J2006h and in lKlvpin et alj i2009l) . We present 
here three sets of fully self-consistent models, all of them with 
a live exponential disk and live dark matter halo with the NFW 
(Navarro et al. 1923) density profile (see Table [T). The live halo 
ensures disk-halo angular momentum exchange, which plays an 
i mportant role in the fo rmation and evolution of bars as discussed 
bv lAthanassoulal ( |2002|) . We simulated barred (B) and unbarred (U) 
models with the aim to compare spiral arms potentially triggered 
by a bar against arms triggered by other mechanisms. Models B 
and U were simulated using the ART code. To generate initial 
conditions of these mode ls we used the Je ans equation moments 
method as introduced bv iHemguislI ( 1 19931) . We have also used a 
multimass method to sample the halo particle distribution, which 
allows us to obtain similar results as using a higher number (N^f /) 
of particles, minimizing t wo body scattering as discussed in 
IValenzuela & Klvpinl jBot) . Barred models have a stellar disk and 
total mass similar to the one obse rved for the MW w ith additional 
initial parameters as proposed by Colin et all ( I2OO6I) (model Kct). 
As discussed in jKlvpin etalj|2002l) . the final properties of this 
model, after rescaling, can reproduce the observed quantities for 
the MW. This rescaling process has not been applied since it does 
not affect the disk kinematic properties analysed in this paper. 
Unbarred models have a smaller disk and a massive and highly 
concentrated halo. Additionally we have used the GADGETS code 
to simulate barred models, labeled G, including a bulge component 
with a Sersic profile (Rt = 1.75 kpc, Mt = 8.57 • 10^ Mq). For 
the bulged model s' initial conditions generat ion we have used the 
code described in lWidrow & Dubinskil ( l2005h . The number of disk 
particles in our models range from one to five million particles (see 

Table [B . 

As expected (e.g. lOstriker & Peebles|[l973l) . the models with a 
relatively dominant disk (B and G) rapidly generate a bar with 
trailing spiral arms while unbarred models (U) do not , at least in 
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Figure 1. Structure of the models. Top: Initial circular velocity of models B (left), G (center) and U (right) , computed from the potential field gradients 
using TIPSY package. We show the total rotation curve (black solid), disk (blue dotted), halo (red dashed) and bulge (green dot-dashed) contributions. Center: 
Density distribution of the models Bl (left), Gl (center) and Ul (right) after 900 Myr of evolution. Bottom: Density distribution of the models B5 (left) and 
U5 (right) after 900 Myr of evolution. The black solid line shows the locus of the spiral amis derived using Fourier analysis (m = 2 for B and G models, 
m = 4 for U). Spiral structure rotates clockwise in all models. 



the first ~ 3Gyrs. In the first case, we have a similar halo and 
disk contribution to the circular velocity inside a radia l expon ential 
length and, as described in IValenzuela & KlvpinI l l2003l) . the 
self-gravity of disk density fluctuations (disk modes) dominates 
making the instability grow and generate a bar (Figure [T] left 
and middle panels). This bar mode induces the formation of a 



bisy mmetric spiral structure app arently connected to the bar ends 
(e.g. iBinnev & Trei^^ llOOSL and references therein). There 
is still controversy on the nature of these bisymmetric spirals 
(resonant coupling, manifolds, or others), but what is evident 
in all our models is that these structure dominates in density. 
For the unbarred model (Figure [T] right) the halo contribution 
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to the potential is higher and the disk modes cannot grow so 
easily. Although in both cases the initial velocity dispersion is low 
(Q = 1.2), the higher halo mass concentration in the unbarred 
model prevents disk from having a dominant bisymmetric mode in 
the first Gyrs of evolution, i.e. other modes grow forming a trailing 
3-4 armed structure. Some other structures with lower density 
appear in the external regions of these simulations (both barred 
and unbarred). In the next section we discuss their imprints on the 
frequency space. As expected for MW like galaxies, we also note 
that rotation curves in our set of simulations are rather flat (see top 
panels of Figure[T}. 

In all our models spiral arm structures are observed for at least 2-3 
Gyrs. We show in Figure |2] the temporal evolution of the spatial 
Fourier modes for the three models with one million particles in 
the disk. We see that dominant spiral modes have a recurrent nature 
with periodicities of less than one galactic rotation. The amplitude 
of the spiral arms in the unbarred case is significantly lower than 
in the barred models. 

In Figure [T] we observe that some parameters of the bar and 
spirals are different from Bl (one million disk particles) to B5 
(five million) models. These differences are arising from the fact 
that the number of particles, spatial and time resolution in B5 
model are much better than in Bl, i.e. we are resolving smaller 
wavelength disk modes in the B5 case than in Bl. These modes 
interact with the disk and then, due to the high non-linearity of 
the system, lead the evolution to a slightly different configuration 
(e.g. bar length and speed). However global quantities like circular 
velocity and density profiles are more robust to such effects. The 
situation is well known and it has been reported i n previous works 
dSellwood & Debattist j[2009l : iKlvpin et alj|2009l) . A convergence 
study of this and other models will be presented in Roca-Fabrega 
et al. (in preparation), however, as will be seen in next sections, 
the main results of our simulations are robust across changes in 
numerical parameters: a bar-arm structure which is dominant in 
density plus external and weaker arms are formed in banned disks, 
while in barless models low amplitude arms are found. In the next 
section we will analyse the kinematics of such spirals structures. 



3 OVERDENSITIES AND ROTATION FREQUENCIES 

We use spatial Fourier analysis azimuthally averaged in 
order to trace the density peak of the spirals and the bar 
jValenzuela & Klvpirj[200 3|), working in cylindrical shells equally 
spaced in galactocentric distance. Figure [T] shows an example of 
how well the spatial Fourier method traces the peak overdensities 
up to the end of the dominant structures (m=2 for barred cases and 
m=3,4 for unbar red ones). W e also u sed a density peak method as 
the one used in iGrand et al. I ( l2012allbl) to test that the results are 
method independent. Once the locus of the spiral is derived, we 
use both, finite differentiation among t hree consecutive snapshots 
and t he classical spectrograms method JSellwood & Athanassoulal 
Il986l) to compute the rotation frequency. 

The advantages of using Spatial Fourier plus Finite Differentiation 
(hereafter SFpFD) are that it allows us to compute the rotation 
frequency of a single known mode m structure and for a single 
time instant. As a consequence, SFpFD is able to show us how 
the struc tures evolve with time. On the co ntrary the spectrogram 
method jSellwood & Athanassoul3 Il98d) needs to be applied 
to a large time interval due to the Nyquist frequency limitation. 
Because of that, the results from spectrograms method may be 



contaminated by recurrent arms sequences with low or negligible 
spiral amplitude due to the transi ent nature of the str uctures. This 
drawback is already discussed bv lGrandetai1 l l2012ah . 
A weakness of using SFpFD, independently of having two or 
more spiral arms, appears when two coexisting structures of the 
same Fourier mode at the same radius are present. In this case, 
the method results in an unique structure placed at the average 
angle. Thus, we have to control these cases to avoid a bias in 
the derived rotation frequency. In contrast spectrograms can find 
discrete rotation frequencies from structures coexisting at the same 
radius without any problem. After carefully weighting the pros and 
cons we state that both methods (SFpFD and spectrograms) are 
complementary. 



4 SPIRAL ARM ROTATION FREQUENCIES 

We perform a first kinematical analysis on the dominant modes 
in density shown in Figure |2] m = 2 for the barred (bisymmetric 
spiral) and m = 4 for the unbarred models (4 armed structure). As 
can be seen in Figure [T] the dominant density structures extend up 
to 10 — 11 kpc. In Figure|3]we present the rotation frequency curves 
for the three models with one million disk particles computed 
using the SFpFD method. It includes all timesteps where we can 
ensure the spiral arm is well formed. We empirically establish that 
a spiral is well formed when the amplitude of the dominant mode 
is above 0.7 times the maximum value of this mode in the range we 
study (see thin black dashed lines in Figure|2}. These density maps 
are constructed from the superposition of all the rotation frequency 
curves at timesteps when the amplitude of the mode is above the 
mentioned threshold. These figures show significant differences 
between barred and unbarred morphologies. Note how in the 
barred models (Bl and Gl, computed using ART and GADGETS, 
respectively) the spiral pattern rotates almost as a rigid body, while 
in unbarred models the rotation frequency curves lie on top of the 
rotation of the disk particles, resulting in a spiral mode corotating 
with the disk. 

A more exhaustive analysis is shown in Figure|4] considering mod- 
els with a higher number of disk particles and a different technique 
for the detection of the spiral structure. Instead of working with 
density maps, here we plot single timestep curves computed using 
SFpFD for models with both one and five million disk particles. 
For the B and G models we avoid the central part of the bar, 
for U the central complex region where the 3-4 armed structure 
converges and in both we also avoid the external regions where the 
number of particles is too low. The red-dashed lines correspond to 
two instants of the one million disk particle models (included in 
the density plot of Figure |3) shown as red dots in Figure |2] More 
important, the blue-dashed lines correspond to models with five 
million disk particles, both barred and unbarred, integrated using 
ART. Furthermore, the blue-dotted line corresponds to the barred 
ART five million model for which the spirals have b een detected 
using th e density peak method (similar to the one in iGrand et"al] 
( l2012allbl) . Note how the conclusions reached in Figure |3] are 
well corroborated. We can see spiral patterns rotating almost 
as a rigid body in all our barred models, both using ART and 
GADGETS, using different techniques for the spatial detection and 
considering a different number of particles. We have also verified 
that simulations with 2 ■ 10" disk particles show the same behaviour 
(not included here). A slightly decreasing rotation frequency with 
radius is observed only for the ART model with one million disk 
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Evolution time (Gyr) 

Figure 2. Amplitude evolution of disk modes. Fourier amplitude for modes 2 (red), 3 (blue) and 4 (green) as a function of time, averaged for radius between 
4 and lOkpc of the barred (Bl, left), the bulge-ban'ed (Gl, center) and the unbarred (Ul, right) models (notice change in vertical scale). Thin black dashed 
lines indicate a threshold in ampUtude used to compute Figure[3]using SFpFD method. The black dots indicate the snapshots for which the kinematic analysis 
is shown in Figure fusing SFpFD method, and thick black dashed hnes indicate the temporal range used to compute spectrograms in Figure |5] 




3456789 10 11 12 3 4 5 6 7 8 9 10 11 3456789 10 11 12 

R (kpc) R (kpc) R (kpc) 

Figure 3. Rotation frequencies as a function of radius calculated using SFpFD method for dominant mode across time in models Bl (left, Gl (center) and Ul 
(right) (see Figure|2). Here we plot the frequency density map of the rotation frequencies computed using all time instants when the spiral arms' amplitude is 
above 70% of the maximum mode amplitude (dashed line in Figure|2). Circular frequencies of disk particles are indicated as sohd white lines and have been 
computed in an intermediate instant of the analyzed time interval. The length of the bar is ~ 4.5kpc and ~ 5kpc for Bl and Gl models respectively. 




3 4 5 6 7 8 9 10 11 12 3 4 5 6 7 8 9 10 11 3 4 5 6 7 8 9 10 11 12 



R (kpc) R (kpc) R(kpc) 

Figure 4. Rotation frequencies as a function of radius calculated using SFpFD method for selected time instants when the amplitude of spiral arms is maximum 
(black dots in Figure|2). All the rotation frequencies analysis of barred (Bl, and B5 left), the bulge-barred (Gl, center) and the unbarred (Ul, right) has been 
computed taking cylindrical shells of O.Skpc width. In red dashed we show results obtained for models with 1 million disk particles. In blue dashed, those of 
the ART models, with 5 million disk particles (left and right panels respectively, for models B5 and U5) when the amplitude of the spirals arms is maximum 
(0.75 and 1.12 Gyr). Circular frequencies of disk particles are indicated as solid lines (red for model s Bl and Ul and blu e models B5 and U5). Additionally, 
blue dots in the left panel show the results of applying a density peak method similar to one used in iGrand etal.N2012ah to B5 simulation, to find the spiral 
structure (a detailed description of this method will be included in Roca-Fabrega et al. in preparation). The dispersion on rotation frequency profiles due to its 
computation at several time steps when the amplitude of the dominant mode and its behaviour can slightly change, is shown at Figure|3] The length of the bar 
is ~ 4.5kpc, ~ y.Okpc and ~ 5.0kpc for Bl, B4 and Gl models respectively. The differences in bar properties between models Bl and B5 are discussed in 
Section |2] 
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Figure 5. Spectrograms for the barred model B 1 (left) and unbarred U 1 ( right) obtained using lSellwood & Athanassoulal il986l) method. We show the Fourier 
component for the m=l,2,3 and 4 modes in a window spanning 0.5 Gyr centered at 1.4 Gyr (left) and 1.45 Gyr (light), (see dashed dark lines in top panels of 
Figure|2j, and with a Nyquist frequency ~ 100 — ISOkms^^kpc^^. The y axis is the angular frequency in kms~^kpc~^ and the x axis is radius in kpc. 
Overplotted in white we show the rotation curve of disk particles computed for an intermediate instant of the analyzed time interval. 



particles (Bl). Even though, this behaviour is completely different 
from all our ART unbarred models, both using one million and five 
million disk particles. 

To ensure that our results are not dependent of the method used 
to compute the rotation frequencies, in Figure |5] we perform one 
last test using, in this case, the classical spectrogram method in 
Bl and Ul models with a Nyquist frequency of ~ 100 — 150 km 
s^^ kpc^'^. The panels framed in red correspond to the dominant 
modes discussed above, and they are to compare with the left 
and right panels of Figure [5] Note that, again, the barred model 
presents a dominant mode that rotates as a rigid body, while in 
the unbarred case, the structures corotate with the disk particles. 
Two features observed in this dominant modes deserve special 
attention. First, in the barred model, a well defined structure is 
observed at a radial range of 11 — 14 kpc (only the beginning 
of the structure is shown in the figure). It rotates slower than the 
disk with a frequency of about 8 km s^^ kpc^^. As mentioned in 



Section 3, here we note one of the advantages of the spectrograms. 
They are able to detect multiple spiral patterns at a given radius. 
Using SFpFD, we have checked that the low frequency structure 
comsponds to one of the lower density structures observed in the 
outer radii of the center-left panel of Figure [T] The complex link 
with the dominant structure will be discussed in a forthcoming 
paper. As it is known, these multiple pattern speeds are common 
in s piral barred galaxies both in simulations a nd observations 
(e.g. iMasset & Tagged [l997l : iButa & ShandllOllI) . Secondly, and 
less important, we can appreciate a flat structure in the rotation 
frequency in the central region of the unbarred model (R < 4 kpc, 
bottom right panel of Figure |5), which could remind the signature 
of a weak bar mode. This is a misinterpretation since it corresponds 
to a m = 4 mode, and, actually, the amplitude of the m — 2 mode 
in this region is less than 2% (see Figure|2j. 

Figure |5] also includes t he spectrograms for the subdominant 
modes. As in Figure 6 of iGrand et al.l (l2012ah we find all these 
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subdominant modes clearly corotating wit h the disk. This be- 
haviour is also observed in other studies as in lMinchev et al. 1 12013) 
(Figure 8), although not discussed there, while they are centered 
on the study of resonant coupling. 



5 DISCUSSION AND CONCLUSIONS 

The rotation frequencies of the spiral modes in barred and un- 
barred models, integrated with different N-body codes (ART and 
GADGET), and analysed with different techniques (spectrograms 
and finite differences), present a well defined different behaviour. 
Whereas unbarred models show transient arms corotating with 
dis k particles, i n good agree ment with t hose recently reported 
by iGrand et al.l ( 1201 2ah and iBaba et all ( 1201 3h . barred models 
deserve further comments. As shown in Figures |3] and |4] barred 
models present a spiral pattern speed almost constant in radius 
for all the range where the spiral structure is dominant in density. 
These results are consist ent with those of models I, II and III in 
iRautiainen & Said ( Il999h . computed using a simplified 2D model 
and with rigid or no halo. We want to emphasize that whereas in 
those early models it is difficult to detect spiral structures in the 
configuration space, particularly in the external regions, in our 
simulations with I and 5 million disk particles, the spiral structure 
is clearly identified (see Figure[T}, allowing us to use the arm phase 
in order to test claimed biases in the derivation of the rotation 
frequencies using spectrograms (see Section|3). 
It i s important to comp are our results with those recently obtained 
by iGrand et al.] ( l20I2bl) . The authors have analyzed N-body/SPH 
simulations of isolated barred galaxies concluding that, spiral 
arms' pattern speed decreases with radius, closely in corotation 
with disk particles. Gran d et a l. (2012b) computed the rotation 
frequencies averaging the values obtained for several snapshots 
over two time intervals during spiral arm evolution, one when the 
bar is well defined in their images (~ 1 Gyr of evolution, see their 
Figure 9) and the other when the bar has significantly weakened 
(~ 1.5 Gyr, Figure 10). In both cases, the spirals pattern speed 
is almost corotating with disk particles. The previous result is at 
odds to what is obtained in our study. However they also noticed 
a small offset of the spiral arms' rotation frequency profile from 
being corotating with disk particles, particularly when the bar 
is stronger, suggesting that further analysis is required. Here we 
have analyzed the correlation between the strength of the bar and 
the spiral rotation frequency in our B I ART model, a Milky Way 
like simulation. With such a purpose, we have computed arms' 
rotation frequency during the first evolutionary stages, when the 
bar is still growing and weak {A2/AQ < 0.1) at about 0.3 Gyr of 
evolution. Only when the bar is weak we observe that the external 
m=2 mode (bisymmetric arm) rotation frequency is close to disk 
corotation. After few galactic rotations, when the bar has fully 
formed (A2/A0 > 0.4) (as shown in Figure |2] left panel), the 
rotation frequency becomes almost constant, approaching the bar 
rigid body rotation, if any with a small decrement in radius (see 
Figure [3] and |4] left and central panels). This result is qual i tative i n 
agreement with the statement discussed by iGrand et al] ( l20I2bl) . 
unfortunately we d o not have information of bar amplitude in 
I Grand et al. I ( l2012bl) simulations, and their models assume a gas 
component making harder a further comparison. We conclude that 
in barred models the spiral rotation frequency approaches the bar 
rigid body rotation with the increment of bar strength. 
Furthermore, our simulations, with a high number of disk particles. 



allow us to analyze the behaviour at larger radii ~ 12 — 14 kpc. 
As can be seen in Figure |5] (m=2 mode, left panel second row), 
the spectrograms method applied to the Bl bar model shows an 
external structure, at radii larger than R= II kpc, that rotates 
slower than the density dominant mode. This external structure 
has a lower amplitude, lower rotation frequency (less than 10 
km s^^ kpc^^), it is tightly wound, and it is recurrently con- 
nected/disconnect ed with the inner and fast er spiral (see movi^B)- 
Other authors (i.e ISellwood&Sparkelll988b have reported similar 
structures but apparently they detected them only in the frequency 
space, probably due to the statistical fluctuati ons given by their 
numbe r of particles (see however figure 16 in IRautiainen & Said 
(Esii)). The lack of detection in configuration space led the 
authors to conclude that inner fast rotating structures correspond 
only to the bar structures while the slower and external ones to the 
entire and unique spiral mode. Here we confirm the detection and 
the rotation frequency in configuration space (see Figure [T] for all 
these structures, and that the inner fast rotating structure includes 
both: a bar and an inner spiral structure apparently connected. The 
slow rotating mode is an external spiral mode. We stress that we 
also find both structures using the SFpFD method in configuration 
space. In conclusion, at large radii, and for m=2, we do see the 
external slow and low amplitude arm coexisting with the dominant 
spiral mode which is most of the time connected to the bar. 
The results presented here for the barred case are also compatibl e 
with the observations recently reported by iMeidt et all (1200 9h . 
They derive rotation frequencies in four spiral galaxies: two strong 
barred, one unbarred and one with rings. They find that in the two 
strong barred cases the inner bisymmetric spirals rotate with the 
same frequency as the bar. Furthermore, one of the banned cases 
present an external armed structure rotating with a much lower 
frequency (see their Figures 8 and 12) as it is observed in our 
BI baiTed model. For the other cases they find multiple pattern 
speeds that can be or not in a resonance coupling situation. It is 
important to mention that at least one of their barred galaxies 
is likely in interaction (MlOl). As discussed by the authors, the 
limiting observational accuracy - precision in the radial binning 
and frequency sampling - do not allow them to confirm or refute 
the small radial decrement in rotation frequency shown in our 
Figure [3] therefore they conclude that arms in their study are rigid 
body rotators. Further observations will be required to confirm this 
tendency. 

From the simulations we have performed it is clear that a difference 
in the spiral arm kinematics exists if they are triggered by a bar 
or by another mechanism, and also that barred models seem to 
show at least two kinds of spiral arms. Further analysis of the arm 
nature, and numerical convergence tests will be discussed in a 
forthcoming paper. 

As a summary we confirm that: 

• The dominant spiral mode (m = 2) in strong barred models is 
most of the time connected to the bar. Its rotation frequency is near 
to the bar solid body rotation, with a small decrement with radius 
which is more important as the bar amplitude decreases. Arms are 
corotating with disk particles only for very weak bars. 

^ This movie can be downloaded from "http://www.am.ub.edu/~J 
sroca/Nbody/movies/Bl.mpeg and it shows the density evolution of 
model Bl, spanning from 0.1 to 3.1 Gyrs. High density regions (~ 3.7 
M0/pc'^) are shown in dark-blue colors while zero density ones are in 
white. 
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• Although the dominant spiral disk mode in barred simulations 
is the one connected with the bar, we observe at least one subdom- 
inant slow and winded mode at large radii, using either SFpFD, 
spectrograms (see Figure |5) and density maps. 

• In unbarred models the spiral structures are cor otating with 
disk m aterial, in agreement with previous results iGrand et aU 
( l2012al) . 

• Our results are robust to changes in numerical parameters 
(time step, spatial resolution, number of particles), mode analy- 
sis techniques (SFpFD, spectrograms, density peak) and also to 
changes in numerical codes (ART and GADGET3). 

A natural question is to ask which is the situation for the MW: 
Do the traditional spiral arms correspond to the mode coupled with 
the bar? or instead correspond to the modes we observed at larger 
radii? or both? Planned and current surveys measuring stellar 
kinematics and distances inside our Galaxy as Gaia (ESA) or 
APOGEE (SDSS) will open up the possibil ity of direct estimation 
throug h methods like the one proposed by iTremaine & Weinberg! 
( Il984l) . Spectro scopic high resolutio n surveys of external galaxies 
like CALIFA l lSanchez et all 1201 2h or MANGA(AS3/SDSSIV) 
will contribute to test our predictions using stellar kinematics. 

We thank A. Klypin, A. Kravtsov and V. Springel for provid- 
ing us the numerical codes and L. M. Widrow for providing the 
code to generate the initial conditions. We thank HPCC project and 
T. Quinn for the implementation of TIPSY package. Finally we also 
thank Luis Aguilar, Ivanio Puerari and Gene G. Byrd for his helpful 
comments on this work. This work was supported by the MINECO 
(Spanish Ministry of Economy) - FEDER through grant AYA2009- 
14648-C02-01, AYA2010-12176-E, AYA20I2-39551-C02-01 and 
CONSOLIDER CSD2007-00050. SR was supported by the MECD 
PhD grant 2009FPU AP-2009-1636. Simulations were carried out 
using Pakal, Abassi2 and Atocatl at lA-UNAM, and Pirineus at 
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